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ON PATHWISE UNIQUENESS FOR REFLECTING BROWNIAN 
MOTION IN C^+'f DOMAINS^ 

By Richard F. Bass and Krzysztof Burdzy 
University of Connecticut and University of Washington 

Pathwise uniqueness holds for the Skorokhod stochastic differen- 
tial equation in C^'^"' domains in M.'' for 7 > 1/2 and d > 3. 

1. Introduction. We will prove pathwise uniqueness for the Skorokhod 
equation in every C^"*""^ domain D C M"^, with 7 > 1/2 and d > 3. We begin 
by giving rigorous definitions of these terms and state our main result. 

We will write x = (xi, . . . , Xci-i,xj) = {x, Xd) for x G M'^. We will say that 
a function <1> iM"^"^ ^ M is C^'^'^ if <I> is bounded and for some constant ci 
and all x and y, 

|V^>(x) - V^>(y)| <ci|x-yP, £,y gM''-^ 

A C^~^^ domain is one which can be represented locally as the region above 
the graph of a C^^^ function. 

Reflecting Brownian motion in a Euclidean domain with Lipschitz bound- 
ary can be represented as a solution to the following Skorokhod stochastic 
differential equation: 

(1.1) Xt = Xo + Wt+ f n{Xs) dLs, 

Jo 

where 11(3;) is the inward pointing unit normal vector at x £ dD, W = 
{Wt - t >0} is a d-dimensional Brownian motion with respect to a filtered 
probability space {Q,J^,{J^t},^), Xq G D, and we require Xt £ D for all t. 
Moreover, Lt is a nondecreasing continuous process that increases only when 
Xt G dD, that is, l^XteD} dLt = 0, a.s. 
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Theorem 1.1. If X and X' are two solutions to the Skorokhod equa- 
tion with respect to the same Brownian motion (hut with possibly different 
filtrations {J-t} and {J^/}, resp.) in a C^"*""^ domain D C M"^, with 7 > 1/2 
and d>3, then 

F{Xt=X[for allt>0) = l. 

Weak existence and weak uniqueness for (1.1) are known; see, for example, 
Theorem 4.2 in [6]. This and the pathwise uniqueness proved in Theorem 1.1 
imply strong existence by a standard argument (see [10], Theorem IX. 1.7). 

Pathwise uniqueness for (1.1) has been proved for all Lipschitz domains 
D cM? with Lipschitz constant 1 in [5] . See [4] for a shorter proof. In par- 
ticular, pathwise uniqueness holds in C^"'"'^ domains in for every 7 > 0. 
See [3] for an introduction to the theory of reflecting Brownian motion and 
[5] for the history of the problem considered in this paper. It was asserted in 
[6] that pathwise uniqueness holds in C^^^ domains for all 7 > 0. The proof 
given in that paper is incorrect; see Remark 3.8 for further details. 

We have a heuristic argument showing that pathwise uniqueness fails in 
some C^'^'^ domains for some 7 > and some d > 3. At this time, we are not 
able to supply all the details needed to turn this claim into a rigorous proof. 
However, we present a counterexample showing that 7 = 1/2 is the critical 
value for Proposition 3.4, a result on pathwise uniqueness for reflected SDEs 
in a half-space that are closely related to the reflecting Brownian motion 
in C^'^'^ domains. The counterexample gives credence to our belief that 
Theorem 1.1 is sharp, that is, it does not hold for 7 < 1/2. 

The proof of Theorem 1.1 is based on an analysis of the distance between 
two solutions to (1.1). We use Ito's formula to express the distance as the 
sum of a martingale and a process with finite variation. The size of the 
oscillations for each of these processes is bounded above using estimates for 
the Green function of the solution to (1.1) and bounds for partial derivatives 
of a mapping of -D to a half-space. Some of the crucial estimates are taken 
from the literature on PDEs in C^"'"'^ domains. We conclude that the distance 
between the two solutions has to be 0, hence pathwise uniqueness holds. 

The paper has three more sections. The next section gives some estimates 
for solutions to partial differential equations. The proof of the main theorem 
is given in the third section. A counterexample showing that Proposition 3.4 
is sharp is given in the last section. 

2. Estimates for solutions to elliptic PDEs. First, we introduce some 
notation. By standard techniques, it is enough to prove uniqueness in the 
region above the graph of a C^'^'^ function Suppose ^iIR'^"^ — > M, and 
D = {{x,Xd):Xd><l>ix)}. 
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Let 

Q{x,r) = {y:\yi -Xi\ < r, i = 1, . . . , d}, 

B{x,r) = {y G M'^ : |y - x| < r}, 

B{x,r) = {y e M'^^^ : |y - x| < r}. 

Let r : ^ M'^ be defined by 

T{x) = {x,Xd-'^ix)). 

We use U = {x^ > 0} for tlie upper half-space. 
Let 

|V$(x) - V$( 



|$|ci+7=sup- ^.^ 
x^y \x-y\^ 

the usual seminorm for the space 

We will use the well-known fact that if h is harmonic in a ball B{x,r) 
and bounded by M in absolute value in B{x,r), then 

(2.1) \Vh{x)\<^; 

see [1], Corollary IL1.4, for example. By induction, the jth-order partial 
derivatives of /i at x are bounded in absolute value by cjM/r^ . 

Proposition 2.1. For any cq <oo there exist ci,C2 < oo such that the 
following is true. Suppose that $(0) = 0, V<^(6) = 0, sup£|#(x)| < 1/4, 
l^lc^+T — CO) cLfid that <l>(x) = if \x\ > 3/4. Suppose h is harmonic and 
bounded in absolute value by 1 in D n Q{0, 1) and h is on each face 
of d{Q{0, l)nD)\ dD. Suppose either (i) /i = on Q(0, 1) n dD, or (ii) 
a/i/9n = a.e. on Q(0, 1) n (9L>. Then: 

(a) we have 

\Vh{x) - Vh{y)\ <ci|x-yP, x, y G 5(0, 5/8) n A 

and \Vh{x)\ < ci /or x G 5(0, 5/8) r\D; 

(b) for any i,j < d, we have 

d'^h 

' " <C2dist(x,9L»)'^-\ X G 5(0,1/2) nZ). 



dxi dxj 



(x) 



Proof, (a) Case (i) follows immediately from Theorem A of [7], if we 
take A{x,y,p) =p and B{x,y,p) = in their theorem. 

Let us consider case (ii). First we claim that h is C°° in Ai UA2, where 

Ai = (5(6, 15/16) \ 5(0, 13/16)) x (0, 7/8), 

^2 = 5(6, 15/16) X (1/2,7/8). 
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Moreover, we claim that the absolute values of the kth partial derivatives 
of h are bounded in Ai U A2 by constants depending only on k and cq but 
otherwise independent of $ and h. 

The result for A2 follows from (2.1) since h is bounded in absolute value 
by 1 and A2 is a positive distance from the boundary of D. We turn to Ai. 
Consider any vertical line segment in Ai of the form L = xq x (0,7/8). Let 
T = i?(xo, 1/32) X (—15/16,15/16) be a tube about this line segment and 
define f{x) on dT to be equal to h{x, \xci\)- If Z is a standard d-dimensional 
Brownian motion and t(T) is the first exit of this process from T, then the 
function 

(2.2) /i(x) =E-7(^.(T)) +E^"'-"'*V(^.(T)) 

is easily seen to have the same boundary values as h on Dn dT. The function 
h is harmonic because each term in (2.2) is separately harmonic in T. We 
have dh/dn = on TOdD by symmetry. Hence h = h inT OU . This implies 
that h is C°° on Ai. The estimate (2.1) and the remark following it can be 
applied to h inside B{x, 1/32) for x € L. We see that the absolute values of 
A:th partial derivatives of h are bounded otl L by constants depending only 
on k. The same holds for h, because h = h on L. This finishes the proof of 
our claim about the behavior of h and its derivatives on Ai yj A2. 

Since h and dh/dn are bounded in Di = (^(0,7/8) x (-7/8, 7/8)) n L>, 
by Green's first identity, 

/ \Vhix)\^dx=\ h(x) — (x)p(dx) < 00, 
JDi JdDi on. 

where p is surface measure on dDi. Hence the function h is in W^''^{Di), the 
Sobolev space of functions whose gradient is in L^. Choose a C°° function 
g : D ^ M of compact support that equals ^ in L>i n -D and has zero normal 
derivative on dU \ Di. Such a g is easy to find because dD is flat outside 
of Di. 

Let / = Ag. We now apply Theorem 5.1 of [9] in the domain {x:xd > 
$(£)}. We set C(s) = sT, so that /(C)(t) = /q C(s)f = cgt^- We set A(x, z,p) = 
p so that diy A{x,u,'Vu) = divVu = Au. We set B{x,z,p) = —f{x), and 
therefore u = g solves div A{x, u, Vn) + B{x, u, Vn) = in the region above 
the graph of We set 'ip{x,u) = 0, so the boundary condition in the the- 
orem in [9] becomes dg/dn = 0. Set = 7. With these choices, verifying 
assumptions (5.1)-(5.6) of Theorem 5.1 in [9] is routine. According to a re- 
mark in the second to last paragraph on page 98 of [9] , one can take f3 = 1 
in Theorem 5.1 of that paper. The conclusion of part (a) of our lemma in 
case (ii) follows now from (5.7) in [9]. 

To prove (b), fix i and let g{x) = ■§^ix)- Then g is harmonic in D and 
is C"> in 5(0,1/2) nL>. Fix xeB{0,l/2)nD and let r = dist{x,dD)/A. Let 
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If he a nonnegative function supported in .6(0, 1) with J ip = 1, and let 
fr{y) = r~''-ip{y/r). Let gr = g*ipr- Since 



\giy) -9r{y)\ 



[giy) - giy - z)]^Priz) dz 



< C3 / Iz]"' ipr{z) dz < C4r^ , 



then by (2.1) 



(2.3) |v(5-g,)(x)|<^ sup \giy)-gr{y)\<cer^~'. 

^ yedB{x,r) 

Note that |V(^r(-2)| < cyr"*^"-^. Using dominated convergence to justify in- 
terchanging differentiation and integration and the fact that / (pr{z) dz = 1, 
we have / \/ipr{z) dz = 0. Therefore 



\Vgrix) 



[g{x - z)- g{x)]Vipr{z) dz 



(2.4) 



< C8 y |2;|'''|V(/3r.(z)| 



Combining (2.3) and (2.4), we obtain 

|Vff(x)|<(c6 + C9)r^-S 

which proves (b). □ 
3. Pathwise uniqueness. 

Lemma 3.1. Suppose € dD and ^>(0) = 0. There exist a continuous 
map H = {hi, . . . , h^) -.D^U with H{0) = 0, e > and ci, C2 > such that: 

(a) H is one-to-one on B{0,e) DD, 

(b) H maps dDnB{0,e) to dU , 

(c) the functions h^, k = 1, . . . ,d, are harmonic on B{0, e) n D, 

(d) hd{x) = OforxedDnB{0,e), 

(e) Vhk -11 = on dDn B{0, e) for k = I, . . . ,d - 1, 

(f ) Vhd{x) = c(x)n(x) for xedDn B{0, e), 

(g) the Jacobian of H is bounded above by 2 and below by 1/2 on B{0,e)n 



(h) dhd/dxd>ci and 

(i) for i,j,k<d, 



(3.1) 



dxi dxj 



(x) 



< C2 dist{x,dD) 



7-1 



xeB{0,e)nD. 
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Proof. It is routine to construct the following functions and the 
corresponding domains D^. For each e > let agree with <I> on B(0,£), 
= if |x| > 3e, and |$£|(;7i+7 < 2|<I>|ci+7. Note that for some 6 > and 
all e € (0,(5), ll^elloo < Let 

= {x: \xi\ < 8e for i = 1, . . . jd — 1 and ^dx) <Xd< 8e}. 

Let be the dilation of by the factor {8e)~^ , that is, = {x :x/{8e) G 
De}. Let 9lZ), = {rr G : 8exd = $.(8ex)}. 

For each e, we define /if, z < d, on Dg, as follows. For z < d, let hf be the 
harmonic function whose boundary values are Xi on dDg; \ dLDi;, and whose 
normal derivative on diDf, is 0. Define /i^ to be the harmonic function whose 
boundary values are x^ on dDg, \ OlD^ and on diDi,. 

By Proposition 2.1, we see that the /if and their gradients are equicon- 
tinuous in i3(0, 1/2) H Dg. Taking a sequence — > along which the 
and their gradients converge, it is clear that in the limit, we obtain for 
i < d a function that is harmonic in Q{0, 1) U, has boundary values Xi 
on (9(5(0, 1) n U, and has zero normal derivative on Q{0, 1) dU. Therefore 
the limit must be the function Xi. For i = d, a, similar argument shows that 
the limit is the function Xd- The gradients also converge, so the limit of the 
gradients of /i^^ must be Cj, the unit vector in the Xi direction. 

Therefore, if we take j sufficiently large, and set e = £j, then the function 

H{x) = {hl'{x/{8ej)), . . .,h'j{x/{8ej))), x G B{0,ej)nD, 

will have all the properties stated in the lemma on B{0,ej) D D. It is routine 
to find an extension H of H to D satisfying the desired properties on all 
of D. □ 

Let E C M*^ be a domain and Z a continuous process taking values in E. 
Let us say Lf is a local time for Z on dE if Lf is a nondecreasing function 
that increases only when Z G dE, that is, /q°° l^ZteE} ^-^f = 

Lemma 3.2. Suppose that H and e > are as in Lemma 3.1, xq G D fl 
-6(0, e/2), and X is a solution to 

Xt = xo + m+ f'n{Xs)dLf, 
Jo 

where W = {W^, . . . , W^) is a d-dimensional Brownian motion and Lf is a 
local time of X on dD. Assume that these processes are defined on a filtered 
probability space {Q.,T,{Tt},^), and W,X and L^ are adapted to {J^t}- Let 
Y = {Y^,...,Y'^) = H{X) and let 

T = inf{t >0:Xt^Dn B{0, 3e/4)}. 



REFLECTING BROWNIAN MOTION 



7 



Then for t <t and i = 1, . . . ,d — 1 




and 








where LJ is a local time ofY on dU that is adapted to {J-t}- 

Proof. The lemma follows from Ito's formula, but we cannot apply 
Ito's formula directly because the function H is not necessarily of class 
on D. For 6>0, let 



and Y^ = {Y^'^,...,Y^^'^) = Hs{X). No te that for 6e {0,e/16), the compo- 
nents of Hs are harmonic and on D n B{0,7e/8). Let = (0, ...,0,1) 
and n{x) = (ni(x), . . . ,nrf(a;)), the inward pointing normal vector on dD. 
By Ito's formula, for t <t and i = 1, . . . ,d, 



By Proposition 2.1, the first partial derivatives of /ij's are Holder continu- 
ous with exponent 7 > 1/2. This and the fact that X is a continuous process 
imply that for any i and k, the processes ^^{Xg + Sea) converge uniformly 

to ^^(Xs) on the time interval [0,r], a.s., when (5 — > 0. It follows that, a.s., 
when 5 — > 0, 



By convention, let n(x) be the zero vector for x ^ dD. Recall from Lemma 
3.1 that V/ifc • n = on 51) n B{0,e) for A; = 1, . . . , d - 1 and Vhd{x) = 
c{x)n{x) for x € dD n B{0,e). We use the Holder continuity of the V/ij's 
and the continuity of X to conclude that for i = 1, . . . , d — 1 and all k, 



Hs{{xi,.. .,Xd)) = H{{xi,.. .,Xd-i,Xd + 5)), 




(3.2) 




(3.3) 



(3.4) 



dhi 



iXs + 6ed)nkiXs)^0, 



dxk 
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uniformly on [0, r] , a.s., when 5 — > 0. This also holds for i = d and d. For 
i = k = d, we obtain 

f)h. 

(3.5) ^(X, + Sea)nk{Xs) ^ c(X,)nfc(X,), 

where c(x) is a bounded function. 

Since H and X are continuous, Yf''' — Yq'^ — > 1^ ~ ^o- We combine this 
observation with (3.2)-(3.5) to see that for t < r and i = 1, . . . ,d — 1, 

(3.6) y/-yj= / Y.^i^s)dw^ 

"'O k=l "-^^ 

and 

(3.7) Y,^ - Y,^ = P^iXs) dW^ + f n{XMXs) dlf . 

Jo ^ OXk Jo 

Note that LJ = Jq c{Xs) dL-^ satisfies the definition of local time of Y on dU, 
that is, it is a nondecreasing continuous process such that I^YsGD} dLj = 
0, a.s., and is adapted to {Tt}. Since H is one-to-one, we can replace Xg 
with H-^{Ys) in (3.6)-(3.7). The last two observations and (3.6)-(3.7) yield 
the formulas stated in the lemma. □ 

Let aij (x) = {dhi / dxj ) o H~^{x) , where the hi 's and H are as in Lemma 3. 1 . 



Lemma 3.3. There exists ci such that for each i,j,k < d, and for each 
xeH{B{0,e/2)r\D), 



(3.8) 



dxk 



< Cl\xd 



17-1 



X € . 



Proof. In view of (3.1), all we have to show is that there exists C2 > 
such that for all x G B{0, e) n D, 

hd{x) > C2 dist(x, (9D). 

But this is immediate from Lemma 3.1(h) and the fact that = on dD fl 
5(0,e). □ 



Let 
(3.9) 



G{x) = \xd 



17-1 
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Proposition 3.4. Suppose a maps to the class of d x d matrices, 
7 > 1/2, and there exists ci such that for each i,j,k < d, 

<Cl\xdV~^, x(^W^. 

Consider the stochastic differential equation 

(3.11) dYt = a{Yt)dWt + m{yt)dL( , ^0 = 2/0, 

where Yt is a d-dimensional process, Wt is a d-dimensional Brownian motion 
with respect to a filtered probability space, yo gU, m is the inward pointing 
normal vector on dU, and LJ is a local time of Y on dU . If Y and Y' 
are two solutions with respect to the same Brownian motion and for some 
bounded stopping time r, 

(3.12) E r{G(Ysf + G(y;)2) ds < oo, 

Jo 

then 

F{Yt = Y/ for allte[0,T]) = l. 

Proof. Let 5 > be small. We can find C2,C3 < oo such that for each 
e > there exists a function : M'^ ^ M such that /g (x) is a nonincreasing 
function of fe{x) = — log |x[ for > e, fe{x) > — loge if |x| < e, and 

(3.13) |V/,(x)|<^, 
\x\ 

We see that df£{x)/dxd < for x e C/ and df£{x)/dxd > for x G C/*^. 
Using (3.10), for any i,j the mean value theorem tells us that 

\aij{x) - CFij{y)\ < C4\x - y\{\xd['~^ V \ydV~^) 

(3.14) 

<Ci\x-y\{G{x) + G{y)). 
Let Zt = Yt — Y/ and apply Ito's formula with the function /g. We have 




(3.10) 



dxk 



(x) 



O'fe 



dxj dx. 



ix) 



< 



C3 



X 



i,j<d, X £ . 
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(3-15) +\fjljP^{Z,) 



d 

Li=l 



ds 



+ / ^{Z,){dLy-dLr) 

Jo OXd 

= MZo) + Mt + At + Vt. 

Note that if G Of/, then Zt € C/^ and so > 0, while if £ dU, 

then Zs & U and so ^^{Zs) < 0. Hence T4 is a nondecreasing process. 
By (3.13) and (3.14), 

\Ar\< r -^\Ys-Y^\\G{Ys)+G{Yi)fds 

Jo \Zsr 

<C6 r{G{Ysf + G{Y^f)ds. 
Jo 

Also by (3.13) and (3.14), 

(M). < r - y,?(G(y.) + G(y;))^ 

Jo \Zs\ 

<cs r{G{Ysf + G{Y:,f)ds. 
Jo 

Therefore, using (3.12), 

E|^,|<C9, E{M)r<Cio. 

So there exists iV > 1, independent of e, such that 



and 



P(K|>iV)<5^<|<« 



Since fe{Zo) > log(l/e) and Vt is nondecreasing, (3.15) and the above esti- 
mates imply that except for an event of probability at most 26, we have 

liminf /e(Z^) > liminf /^(Zq) - 2iV = oo. 

We conclude that Z^ = 0, with probability greater than or equal to 1 — 25. 
Since 6 is arbitrary, Zj- = 0, a.s. This is true, a.s., simultaneously for all 
stopping times t At, for rational t. Since Zt is continuous in t, our result 
follows. □ 
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Corollary 3.5. Suppose that e > is as in Lemma 3.1, xq ^ D r\ 
B(0,e/2), and X and X' are solutions to 

Xt = xo + Wt+ f'n{Xs)dLf, 
Jo 

X[ = xo + Wt+ f n{Xi) dLf , 
Jo 

relative to the same Brownian motion W . Let 

T = mf{t >Q:XtiDr\ B{0, 3e/4) or X[ ^ D n B{0, 3e/4)}. 

Then 

F{Xt = X[ for all t € [0, r]) = 1. 

Proof. Let H be as in Lemma 3.1 and Y = H{X) and Y' = H{X'). It 
will suffice to show that P(lt = Yl for all t € [0, r]) = 1. In view of Lemma 3.2 
and Proposition 3.4, all we have to show is that (3.12) holds for these pro- 
cesses Y and Y' . 

Extend the definition of each aij by setting aij{x, x^) = crij{x, for x G 
W^\U. By [8] (see also [2]), the Green function K{x,y) for the corresponding 
elliptic operator is comparable to \x — y\'^~'^. Let Ui = H{B{0,e) (ID). If 
Yo = yo, then 

¥.rG{Ysfds<[ K{yo,x)G{xfdx. 
Jo JUi 

We will estimate the last quantity in the special case yo = 0. The general 
case can be dealt with in an analogous way. Let A^j = {x : < < 
2-fc+i^2-J <Xd< 2-J+i}. For x e A^j, we have K{yo,x) < ci2~''(2-d) ^^^^ 
G{xf < C22-2{7-i)j. The volume of Akj is bounded by 032''''^'^-'^^-^. Since 
7 > 1/2, the series ^j>i2 •'^ ^ 'is summable, and we have 




<C4 5^2-'=^2-J'(27-l)<oo. 
k>l j>l 

It follows that E £ G(y,)2 ds< 00 and similarly, E £ G{Y!.Y ds< 00. □ 

Proof of Theorem 1.1. By standard arguments, the local version of 
pathwise uniqueness proved in Corollary 3.5 can be extended to the global 
assertion stated in Theorem 1.1. □ 

Remark 3.6. The same proof as that in Proposition 3.4 proves the 
following theorem. 
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Theorem 3.7. Suppose Z and Z' are two solutions taking values in 
to the stochastic differential equation 

(3.16) dZt = aij{Zt)dWt, 

where Wt is a d- dimensional Brownian motion. (We allow Z, Z' to he adapted 
to possibly different filtrations.) Suppose there exists a function G:W^ ^ 
[0, oo) such that for each t < oo, a.s., 

l\G{Zsf + G{Z'^f)ds<oo, 
Jo 

and for each i,j = l,...,d and y € M.'^, we have 

(3.17) \aij{x) - a^j{y)\ <\x- y\{G{x) + G{y)). 
Then 

F{Zt = Z[ for allt>0) = l. 

Remark 3.8. The proof in [6] that weak uniqueness holds for (1.1) is 
correct, but the proof of strong existence is not. If {J-t} is the filtration of 
the Brownian motion, it was proved that there exists a solution Xt of (1.1) 
with Xt being J^i measurable for all t <1, but the process X constructed 
there was not necessarily adapted, that is, it was not shown that Xt was J^t 
measurable. 



4. Counterexample. Let U be the upper half-plane and consider the 
stochastic differential equation 

(4.1) dYt = a{Yt)dWt + n{Yt)dLl, Yo = yo, 

where is a three-dimensional process, Wt is a three-dimensional Brownian 
motion with respect to a filtered probability space, yo €U, n is the inward 
pointing normal vector and LY is the local time of Y on dU. 

Theorem 4.1. For every 7G (0,1/2) there exists a which maps U to 
the class of 3x3 matrices such that a is bounded, uniformly positive definite, 
there exists c\ such that for each k <3, 

'7—1 

< ClXg , X GU, 

and there exist two solutions Y and Y' to (4.1) with respect to the same 
Brownian motion and with yo = such that for some t > 0, ¥{Yt 7^ Y/) > 0. 

Proof. We will identify dU and M^. Choose a function : ^ [0, 00) 
with the following properties: 



(4.2) 



dxk 



{x) 
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(i) ^{x) = ipdx]) ior some ip :M ^ [0,00). 

(ii) ^{x) is C°°. 

(iii) ip{0) = 1, V' is decreasing on [0, cxo), ip{r) 



for r >3/4. 



For integer k, let 

(4.3) ^l{x) 

(4.4) ^k{x) 



E 

(ii,j2)ez2 



(3il2-^3f22-*^)). 



Let ifk be the harmonic function in U which has boundary values ^^(2;) 
on dU. Consider a large integer ni whose value will be specified later and 

let (f = J2m.>0 ^mni ■ 

Let aij = for i 7^ j, 0-33 = 1 and (Tjj{x) = 1 + ^{x) for j = 1, 2. 
First we will show that the CTjj's satisfy (4.2). For some 02,03 < 00, all j, 

are harmonic, we have |v9o(x)| < C2 and ||^(x)| < C3 for all x € By (2.1), 
|^^(x)| < C4/X3 for X3 > 1. We use scaling to see that 



and all x £ dU , we have 19^0(2^) I ^ C2 and |§^(a;)| < C3. Since (^o and 



dxj 



x) 



< 



Hence, for some c-j < 00, j 

E -52 



('c52-^'(t-i), ifx3<2-*^, 
\c62-'=V^3, ifx3>2-'=. 
1, 2, 3 and all x E ^7, 



< 



-fc(7-l) 



+ 



k:X3<2- 



E 

fc:X3>2-'= 



7-1 



This and the definition of the Cij's imply that (4.2) holds. 

Assume that pathwise uniqueness holds for (4.1). We will show that this 
assumption leads to a contradiction. Pathwise uniqueness together with 
weak existence imply strong existence. Let Y and Y' be the unique solu- 
tions to (4.1) with the same driving Brownian motion W, with starting 
points Yq and Yq. Most of the time we will be concerned with Y and Y' 
starting from different points, that is, Yq ^ Yq. 

We will write Yg = {Yi{s) ,Y2{s) ,Y3{s)) , and similarly for Y^. 

Note that the third component of y is a one-dimensional reflecting Brow- 
nian motion starting from and the same is true for Y' . Moreover, if 
^3(0) = ^'(0)> then the third components of Y and Y' are equal to each 
other for all t, a.s. 

By Ito's formula, for A; = 1, 2, 

{Y,it)-Yl{t)f 

= (n(0) - Yl{0)f + r 2{Yk{s) - Yl{s)){akk{Ys) - akkiY^)) dW^ 
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+ i f 2{akk{Ys) - akk{Y^) f ds 
Jo 

= (y,(o) - Ymf + r 2(n(s) - - ^(nO) 

/o 

It follows that if we set Rt = \Yt — Y^ , then 

2 







+ / i^{Ys)-^iY:)fds. 



Rt = Ro+ i2'^i^k{s)-Yas)){ipiY,)-ip{Y:))dW^ 
(4.5) 

+ fY.iv{Ys)-^{Y^)fds. 
^0 fc=l 



This implies that Rt is a time-change of the square of a two-dimensional 
Bessel process. Hence, the process t ^ \ Yt — Y/\ has the same exit probabili- 
ties from an interval as a two-dimensional Bessel process. The main technical 
goal of this proof will be to show that R hits 1 before some time to < oo 
with probability po > 0, where to and po are independent of i?o > 0. 

We will now derive some estimates based on (4.5) that will be needed 
later in the proof. 

Suppose that y(0) / Y'{0), YsiO) = Y^{0) and consider 6 > 0, oq € [0, 1) 
and ai > 1. Let 

n = inf{t > 0: l^t - y/l ^ [ao\Yo - Y^lai\Yo - Y^\]}, 

T2 = inf |t > : j\^{Ys) - ^{Y^)f ds > b\Yo - Fo?}. 

Since |y — y'| is a time-change of a two-dimensional Bessel process, we have 
for a two-dimensional Brownian motion W, 



Yo - Y^ + w(^j\^iYs) - ip{Y:)f ds^ 



lYt-Yl 

Note that if T2 < oo, then 

Hi^iYs) - ip{Y:)f ds = b\Yo - Yo?, 
Jo 

and therefore 

\Yr, - y;i = 1^0 - Y^+W{b\Yo - yo?)|. 

Hence, if we take a = (Yq - >o)/l^o -Y^\, 

P(t2 < ri) < P(t2 < oo, \Yr, -Y;je [ao\Yo - Y^\,ai\Yo - Y^\]) 

< F{\Yo - Y^+W{b\Yo - Fof )l G Nl^o - Y^\,ai\Yo - Y^\ 
= P(|a + W(6)|G[ao,ai]). 
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The last quantity is a constant po > depending on ao,ai and b but not 
depending on Y or Y' . Moreover, it is easy to see that for fixed uq and ai, 
Po^O when 6 — oo. We record the above inequahty for future reference as 



(4.6) 



(t2 < Tl) <P0- 



We will now estimate the probability of the complementary event in the 
same setting. We have 

lP(n<r2)<p( U {\Yt-Y^\(^[ao\Yo-Y^\,ai\Yo-Y^\]}\ 

\0<t<T2 I 

<P( U {|yo-Fo' + W^(s|^o-l?,?)|^[ao|>o-l?,'|,ai|>o-^ol]} 

\0<s<6 J 



P U {|a + VF(s)|^[ao,ai]} 

\0<s<6 J 



This easily implies that 
(4.7) 



P(n<r2)<pi, 



where pi ^ as 6 — > 0, for any fixed ag and a\. 

Take 1 — oo = ai — 1>0 and consider an arbitrary finite stopping time r. 
Note that 



|||y, - y;i - 1^0 - y^W > («i - 1)1^0 - Y^\, 



((^(n) - ^{Y:)f ds < b\Yo - ^0? \c{n< T2}. 



i\2 



If we take 02 = ai — 1, it follows from (4.7) that for any 02,^*0 > there 
exists 6 > such that for every stopping time r, 



Yr - Y;\ - \Yo - Y^W > 021^0 - 

(ifiYs) - V9(y;))' ds < b\Yo - Fo? ) < PO. 



(4.8) 



We will next estimate <f{x) — <p{y) for x and y such that X3 = 1/3 > 0. 
The function ipQ is not constant on horizontal planes, so for some £U 
we have (|^, ^){x^) = v / 0. Let ro > be such that S(x°,ro) C U and 



(4.9) 



\dx 



d(po\, . 



< IvI/lOO 
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for all X € B{x^,rQ). By scaling, 



(4.10) 



<2-fc(7-i)|v|/l00 



oxi 0x2 / 

for A; > 1 and all x e B{2-''x^ 

We have [(/^^(x)! < 2~^'^ for all x, so X]fc>o ll'/'fclloo < co- It follows that 

00 

(4.11) lim sup Vvmn(x) = 0. 
By (2.1), 

00 

lim sup V |Vv9m„(x)| = 0. 

"^°°xeS(xO,ro)^=l 

We choose ni so large that 

00 

(4.12) sup |Vc^^n,(x)| < |v|/100. 

xGB(a;0,ro)m=l 

By scaling, 

00 

(4.13) sup |V^.„„i(x)|<2-'=(^-i)|v|/100. 

xGB{2-fexO,2-fero)m=fc+l 

For some c% < 00, by scaling, ||V(^fc||oo < Cs2~''^"'~^\ We make ni larger, 
if necessary, so that 

-1 -1 

(4.14) l|V^„^nJ|oo< E C82— "1(^-1) < |v|/100. 
rra=— 00 m=— cxD 

We use scaling again to see that 

E l|Vv^™„J|oo< E C82— "1(^-1) <2-'=(^-iVl/100. 

m=— 00 m=—oo 

Then 

fc-i fc-i 

(4.15) El|V<^-"ilU^ E l|V¥'mnil|oo<2-'=(^-^)|v!/100. 

m=0 m=— 00 

With this choice of ni, we have by (4.10), (4.13) and (4.15), for all x G 

B{2-''x'^,2~''ro), 



(4.16) < 



(I?' I?) ^ |v^..,(.)|+Ellv^^ 



m=k+l m=0 



^2^2-fc(7-i)|v|/i00. 
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By invariance under rotations by the angle '7r/2, there is a vector w, 
perpendicular to v and to the vertical axis, such that (4.16) holds with 
replaced by a different reference point, say, x^; that is, for A; > 1 and 



X € B{2~'^x 
(4.17) 



,2-«ro) 







\dxi 


8X2/ 



(x) 



W 



< 3. 2-fc(7-i) 1^1/100. 



Estimates (4.10), (4.13) and (4.15) are invariant under shifts by vectors of 
the form (3 • , 0) with (11,12) €Z , so we have for (^1,^2) , 



(4.18) 





dip_\ 


\dxi 


8x2) 



(x) 



'fc(7-l). 



<3.2-fc(7-i)|v|/l00, 



for xG5(2-''x° + (3-2-*^n,3-2-'=Z2,0),2-'^>o), and 



(4.19) 







\dxi 


8X2/ 



(x) 



-fc(7-l) 



W 



<3-2 



-fc(7-l) 



w\/100, 



for X G B{2~^x^ + (3 • 2-*^n,3 • 2-'=i2, 0), 2-*^ro). 

Fix some such that 2~^° < ro/4 and consider any vector u with zero 
third component such that 2~^°~'^~^ < |u| < 2~^°~^ . There exists cg < 00 
such that for every z € 8U we can find x^ and x^ of the form 



2-'=x° + (3-2-''ii,3-2-'=i2,0), 



(4.20) 

= 2-'=x^ + (3 • 2-''n, 3 • 2"''i2, 0), 

with the property that \z — x^| < cg2~'^ and \z — x^\ < c<j2~^ . One of the 
vectors u or — u must form an angle smaller than 37r/8 with one of the 
vectors v or w. Suppose that u forms an angle smaller than 37r/8 with v and 
consider x € i?(x^, 2~''>o/100) and y S B{x'^ + u, 2~''>o/100). By integrating 
Vy5 along the line segment joining x and y, and using (4.18), we obtain 

^{y)\>ciQ2~^\ 



(4.21) 



|V9(x) 



where cio does not depend on k,z,x or y. In the remaining cases, either 
the above estimate holds as stated or it holds for x € i?(x^, 2~'^ro/100) and 
y GB(x3 + u,2-'^>o/100). 

Next we will estimate the rate of growth of t ^ Jq{ip{Ys) — ^p{Y!.)Y ds. 
Suppose that ^0,^0 € dU and \Yq - Y^\ G [2-^-^,2'^]. Let k = j- ko, n = 
ro/lOO and u = Yq —Yq. Assume without loss of generality that u forms an 
angle smaller than 37r/8 with v. Let x^ be defined as in (4.20), relative to 
z = Yq. Then (4.21) holds for x e B{x^,2~^ri) and y G i?(x^ + u, 2~^ri). Let 
r = inf{t > : LJ > 2~^}. Let Ai be the following event: 

The process Y3{t) reaches level 2~^ at some time Ti < r, then Y hits 
B{x^,ri2~^~^) at a time T2 < r, and then stays in i?(x^, 3ri2~^/4) 
until at least r2 + 2"^^. 
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The probability of Ai is bounded below by pi > 0, independent of k. To 
see this, we argue as follows. Extend a to all of by reflection, that is, 
define a{xi,X2,X3) = a{xi,X2,—X3) = a{xi,X2,xs) for rc3 > 0. Let be a 
weak solution to 

(4.22) dVt = a{Vt)dWt, ^0 = (0,0,0). 

The solution to (4.22) is weakly unique since a is continuous and uniformly 
positive definite (see [3], Chapter VI). Let Vt = (yj\y/,y/) be defined by 
Vf = VI for i = 1, 2 and Vt = l^'^l- A routine calculation with Ito's formula 
shows that 

dVt = cj{Vt)dWt + ^{Vt)dLY 

for some three-dimensional Brownian motion W . The proof of the theorem 
of Yamada-Watanabe (see [10], Section IX. 1) applies in the present context, 
so since we are assuming that there is pathwise uniqueness for (4.1), then 
there is weak uniqueness for (4.1). Hence V and Y have the same law, and 
to estimate the probability of j4i, it suffices to estimate the probability of 
Ai when Y is replaced by V . By the definition of y, it is enough to estimate 
this probability when V is replaced by V . This latter probability is bounded 
below, using scaling and the support theorem for diffusions corresponding 
to uniformly elliptic operators; see [3], Theorem V.2.5. This completes the 
proof that P(^i) >pi. 

Consider some positive constants cn and ci2 such that cn < ri/8 and 
cii2-^ < \Yq - yo'l/2, and let 



A2 = {\\Yr- Y;\ - \Yo - Y^W > cii2" 

riipiY,)-^iY^)fds<c,22^ 
Jo 



■2k 



By (4.8), we can choose ci2 > small so that P(^2) < Pi/4. Then P(^i \ 
A2) >Pi/2. Suppose that the event Ai \A2 occurred. Then one possibility 
is that 

(4.23) [{^{Ys) - H^{Y;,)f ds > ci22~^\ 

Jo 

If the event in (4.23) does not occur, then Y^ G B{x'^ + u, ri2"'^) for t G 
[r2,r2 + 2-2'=] and so by (4.21), 

(v9(y,) - (^(y;))2 ds > {^{Ys) - ip{Y:)f ds 

- Jt2 
(4.24) 

>Ci32-2'=^2-2'= = Ci32-2'=(l+^). 
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Combining the two cases (4.23)-(4.24), 

(4.25) P^"((/.(n)-v9(y;))2ds>ci42-2fc(i+7)^ >p^/2. 

Recall our assumption that \Yq — Yq\ G [2^-'^^, 2^-']. Let 

n = mf{t > : - y/l ^ [\Yo - Y^\/2, 2\Yo - Y^\]}, 

n = inf |t > : ^ {ip{Ys) - ip{Y^)f ds > 0152-2^' |, 

for some constant C15. According to (4.6), we can choose C15 sufficiently large 
so that 

(4.26) P(r3>r4)<l/2. 

Recall that k = j — for a suitable constant /cq. Let tq = 0, ri = r A T3, 
Tj = T3 A (r o 6't-,_i + ri_i) for i > 2 and 

= { //'^' (^(^«) - '^(^«'))' ^ ci42-2(^+i)(i+^) I u {t,+i = Ts}, 

for i > 0. The events Fi are not necessarily independent, but if we take 
Ti = (j{Fi, . . . , Fj_i), then we have 

P(i^d-^*-i)>Pi/2, 

by (4.25) and the strong Markov property applied at the r^'s. Hence, we 
can compare the sequence {-Fi} to Bernoulli trials to reach the following 
conclusion. Let /q be the smallest integer such that X]i=o '^Fi — where 
C16 is chosen so that ci622*^'^ci32-2(^+i)(i+^) > C152-2J. Then E/q < 01722^=-^. 
It follows from the definition of Fi 's and constant c\q that r/,, > r4 A r3 . 
Note that L^^, - FX < 2~K Hence, 

(4.27) EL^^/\r3 < EL^^^^ < 2"%/o < 2--'ci722'=^ = cis2~^^^~'^'<\ 

Let 5i = r4 A Ta, 5, = T3 A {T^ o 63^., + Si-i) and h = inf{i : 5, > r3}. 
Let Gi = {T3 < Sj} and Qi = cr{Gi, . . . , Gi). By the strong Markov property 
applied at the stopping times Si and (4.26) we have F{Gi \ Qi-i) > 1/2 for 
all i. Thus we can compare the sequence {Gi} with a sequence of Bernoulli 
trials to reach the conclusion that E/i < cig < 00. 

By (4.27) and the strong Markov property applied at the Si^s, E(L|i.^^ — 

-^5.) < cis2~^^^~'^"'^ for all i. A standard argument based on the strong 
Markov property shows that 

EL^3 < EL^^ < ci82"'^'(^~2^)E/i 

(4.28) 

<ci82-'=(^-2^)ci9 = C2o2-'=(^-2^). 
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Assume that IIq — ^ol ~ ^ some (large) integer j, and 13(0) = 5^3(0). 
Let C/i = Ts, [/, = Ts o eu,_^ + C/»-i for i > 2, and h = mf{i : \Yu^ -Y{j^\> 
1}. Let NP be the number of downcrossings of the interval [2"*^"^, 2"^^], 
that is, is the number of distinct i < I2 such that \Yu^ — YIj. \ = 
and jyi7j_|_i — \ = 2~^~^. Similarly, let be the number of distinct 
i < h such that \Yui - 'i^^l = 2"'="^ and \Yu,+^ - ^L^i+il = 2"'"- The process 
t ^\Yt — Y{\ is a time-change of a two-dimensional Bessel process. For any 
A; > 1, the two-dimensional Bessel process starting at hits 1 before hitting 
2~^~^ with probability 1/(A; + 1), so Nj? is stochastically majorized by a 
random variable with a geometric distribution with parameter + 1). 
Therefore, EAf < C2ik. Since \N]^ - N^\<1, we also have "KNV < C22A:. 
Let Nk be the number of distinct i < I2 such that \Yu^ — Y^. \ = 2~ . We see 
that EiVfc < ENj/^-^ + EiVf < C23A;. 

If lYjy. — y^. I = 2~'^, then by (4.28) and the strong Markov property, 

E(L^,^^ — LJjJ < C2o'^~^^^~'^"'^ ■ It is routine to combine this with the es- 
timate EAfc < C23A; to see that 

Thus EL^^^ < C24 Efc>i A:2~*^(i-27) < (jug ^ the assumption 7 G (0, 1/2). 

Let Ts = inf{t > : |Ft - y/| = 1} and note that T5 = Ui^ . The distribution of 
{LJ ,t > 0} does not depend on {Yq — Yq \ and we have = 00, a.s., so there 
exist P4 > and to < independent of lYo — ^ol' such that P(C//2 < to) = 
P(T5 < to) > p4. In other words, there exist P4 > and to < 00, independent 
of \Yo - Fo'l, such that P(suptg[o,t„] \Yt - y/| > 1) > p^. 

We finish by using the following standard argument. We have assumed 
that there exists a unique strong solution to stochastic differential equation 
(4.1). Let Y be the solution with Yq = (0, 0, 0) and let Y^ be the solution with 
Yq = (0,2~-^, 0). It is easy to see that the family {Y,Y^ ,W}^j>ij is tight in 
the space C([0, 00), M^)^ equipped with the topology of uniform convergence 
on compact sets. Take a weakly convergent subsequence (as j 00) and let 
{Y, Y' , W} be the weak limit. Note that (Y^)^ = Y^ for each j, and therefore 

[\{Yi)dLf= fniY,)dLj 
Jo Jo 

for each j and each t. With this observation it is clear that both Y and Y' 
are solutions to (4.1) with the same driving Brownian motion W and both 
starting at (0,0,0). By the assumption of pathwise uniqueness, If = 1^' for 
all t>Q, a.s. However, we showed that for some tQ,po, 6 > 0, 

p( sup \Yt-Yt^>6) >po 
Vte[o,to] / 



REFLECTING BROWNIAN MOTION 



21 



for all J > 1. Passing to the limit, we see that with probability at least 
Po, we have supjg[o,to] 1-^* ~ \ > ^f^i ^ contradiction. Therefore pathwise 
uniqueness cannot hold for (4.1). □ 
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